Spin-orbit anisotropy measured using ballistic spin resonance 
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Spin relaxation can be greatly enhanced in narrow channels of two-dimensional electron gas due 
to ballistic spin resonance, which is mediated by spin-orbit interaction for trajectories that bounce 
rapidly between channel walls. The channel orientation determines which momenta affect the re- 
laxation process, so comparing relaxation for two orientations provides a direct determination of 
spin-orbit anisotropy. Electrical measurements of pure spin currents are shown to reveal an order of 
magnitude stronger relaxation for channels fabricated along the [110] crystal axis in a GaAs electron 
gas compared to [110] channels, believed to result from interference between structural and bulk 
inversion asymmetries. 



Spin-orbit interaction (SOI) is one of the most promis- 
ing tools for fast spin control in solid state systems pQ. 
Although SOI was first studied in the context of spin 
relaxation and decoherence, a variety of recent propos- 
als and experiments have explored the use of SOI for 
intentional spin separation, filtering and coherent con- 
trol in spintronic devices and spin qubits [2HZ]- Two- 
dimensional electron gases (2DEGs) make an attrac- 
tive platform for spin-based electronics because elec- 
tron transport can be easily controlled using electrostatic 
gates. 

For 2DEGs in zinc-blende semiconductors such as 
GaAs, the primary contributions to SOI are structural in- 
version asymmetry (SIA) and bulk inversion asymmetry 
(BIA), which add to induce an effective magnetic field, 
B so , that acts on the electron spin [HE?]: 

B s ° = —((«- p)k m i ~(a+ R)k lw 3) (1) 

to first order in the momentum k, where i and j are 
unit vectors along the [110] and [110] crystal axes, and a 
and (3 represent the strength of the first-order SIA and 
BIA terms, respectively. SIA and BIA vary from 2DEG 
to 2DEG depending on sample details, but a and j3 are 
believed to be of the same order of magnitude in typi- 
cal GaAs heterostructures (triangular wells). By Eq. 1, 
this should lead to anisotropy in the spin-orbit field, as 
the effects of SIA and BIA add to or subtract from one 
another depending on the direction of the electron mo- 
mentum [see Fig. [lja)]. When a and (3 have opposite 
signs, the interference between SIA and BIA gives rise to 
B so that is stronger for momenta in the [110] direction 
[Fig. [ljb)]; this anisotropy is reversed if a and /? have 
the same sign. 

Measuring spin-orbit anisotropy is crucial to the design 
of electronic devices for spin manipulation in a 2DEG, 
and for optimizing heterostructure growth for spintron- 
ics. The condition \a\ = \/3\ is of great interest for appli- 
cations where one seeks to rotate electron spin using SOI 
with d.c. voltages [TOHTi] . Conversely, spintronics pro- 



posals based on the mesoscopic spin Hall effect depend 
on moving to \a\ <C \f3\ or |a| 3> |/3|[T5], Despite the 
importance of spin-orbit anisotropy to 2DEG spintron- 
ics, there are only a few published reports of a and f3 
in the triangular well heterostructures typically used for 
electronic devices, and the values found in these reports 
are often inconsistent [TMTS] . The only transport mea- 
surement of spin-orbit parameters to date was based on 
indirect technique in which the shapes of magnetoresis- 
tance curves were compared to a detailed theory of weak 
antilocalization[18 . 

In this Letter, we report an all-electrical tech- 
nique for directly accessing spin-orbit anisotropy in a 
GaAs/AlGaAs 2DEG. The method is based on ballistic 
spin resonance (BSR), which suppresses spin relaxation 
length in narrow conducting channels [19] . Spin-orbit 
anisotropy is measured by comparing BSR in 1 /zm-wide 
2DEG channels fabricated along the [110] and [110] crys- 
tal axes. The spin relaxation length at resonance is 3 
fjja. in [110] channels, compared to 40 /<m in [110] chan- 
nels. Based on Monte-Carlo simulations of semiclassical 
spin dynamics, this anisotropy in relaxation length im- 
plies a spin-orbit anisotropy \a — a + f3\ — 15±5, the 
largest value reported so far for GaAs triangular wells 
|20j . Such a strong anisotropy will play an important 
role in the design of spin qubits and spin transistors in 
similar materials. 

The channels are defined by electrostatic gates on the 
surface of a [001] GaAs/AlGaAs heterostructure. The 
electron gas (n s = 1.11 x 10 11 cm~ 2 and mobility // = 
4.44 x 10 6 cm 2 /Vs at 1.5K) is 75 nm below a layer of Si- 
doped AlGaAs, and 110 nm below the wafer surface. The 
data presented in this paper are from four devices: short 
and long channels (length 100 and 130 /im) aligned along 
[110] and [110]. In short channels, the spacing between 
injector and detector contacts is Xid — 5 /im; in long 
channels, Xid = 25 /itm. Measurements are performed 
with an a.c. lock-in technique in a dilution refrigerator, 
with an external magnetic field B ext applied in the plane 
of the electron gas. 
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FIG. 1: (a, b) Spin-orbit fields due to SIA and BIA and the 
total spin-orbit field, for the case a = —1.5/3. (c) Optical im- 
age of a typical device (short channel). In the measurement, 
S-gate is depleted to form two 1 /im-wide channels from this 
single device, (d) Schematic of measurement setup. The 
channel and the injector and detector QPCs are defined by 
depleting the 2DEG with gates (light gray). A-gate changes 
the distance between the equilibrium reservoir and the injec- 
tor, (e) Periodic trajectories of electrons (dashed line and 
thin arrows) between channel boundaries are accompanied by 
an oscillating B so (thick black arrows). Injected spins are 
aligned transverse to the channel by the in-plane field B ext . 
(f ) Component of spin-orbit field along the channel axis os- 
cillates with different magnitudes for [110] and [110] channels. 



Two quantum point contacts (QPCs) are embed- 
ded into each channel. At high magnetic field the 
QPC conductance traces show spin-polarized conduc- 
tance plateaus at 1 e 2 /h, which are used to generate and 
detect electron spins [see Figs, gc) and [jjd)] [2H23]- 
The injector QPC drives polarized electrons towards the 
left end of the channel. Diffusion of spin polarization to- 
wards the equilibrium reservoir on the right end creates a 
pure spin current. Spin polarization accumulated above 
the detector QPC is then quantified by the non-local spin 
signal V nl . 

Electron trajectories in these devices bounce back and 
forth many times between the channel walls before scat- 
tering, because the mean free path is much longer than 
the channel width. Ballistic spin resonance is driven by 



this motion, and mediated by spin-orbit coupling. As 
shown in Fig.[TJe), an injected electron bouncing between 
the channel walls feels a spin-orbit field B so , and the com- 
ponent of B so along the channel axis oscillates with the 
bouncing frequency [Eq. [I] . 

The magnetic field B ext sets the spin quantization di- 
rection and the polarization of injected spins but does not 
affect the orbital motion of the electrons because it is ap- 
plied in the plane of the electron gas. BSR is observed 
when B ext is applied perpendicular to the channel axis. 
In this case, the periodic component of B so is transverse 
to the external field, and it induces spin rotations when 
the Larmor precession frequency g/iB ext /h matches the 
typical bouncing frequency vf/^w, where Vf is the Fermi 
velocity and w is the width of the channel. Spin polariza- 
tion is rapidly lost on resonance because every electron 
follows a different bouncing trajectory. 

Figure 2 shows characteristic traces of V n i(B ext ) for 
four different channels. Two qualitative features are ob- 
served in all data sets. First, V n i starts near zero and 
generally grows with magnetic field, showing the rise in 
QPC polarization as the Zeeman energy becomes larger 
than thermal and tunnel broadening [23] ■ Second, a dip 
is observed between 5 and 8 T in each device, resulting 
from faster spin relaxation at the BSR condition as is 
described below. The centers of the BSR dips in the four 
panels of Fig. [2jvary due to changes in the 2DEG density 
(and hence vp) for different cooldowns. 

Measurement of BSR in two orthogonal channels pro- 
vides direct access to spin-orbit anisotropy. The channel 
orientation determines whether [110] or [110] component 
of momentum oscillates: [110] -oriented channels induce 
oscillating fcj 10 , while [110]-oriented channels induce os- 
cillating fcno. According to Eq. [I] an oscillating fcj 10 
generates an oscillating Bfi Q oc (a — j3), whereas fcno 
generates B{° w oc (a + 0) [Fig. [iff)]. 

The spin signal vanishes almost completely inside the 
dip in [110] channels, while only a weak suppression is 
observed in [110] channels [Fig. [2]. At a qualitative level, 
this indicates that the oscillating spin-orbit field in [110] 
channels is much larger than in [110] channels, and based 
on the direction of anisotropy we conclude that a and 
(3 as defined by Eq. [I] have opposite signs. BSR was 
measured in eight [110] and three [110] channels in this 
2DEG. Details of the V n i(B ext ) traces were different for 
each channel (even each cooldown) due to variations in 
QPC polarization, electron density, etc., but the dra- 
matic difference in relative BSR dip depth for [110] versus 
[110] channels was consistent for all devices. 

A quantitative determination of spin-orbit anisotropy 
is made by extracting spin relaxation length on reso- 
nance for the two channel orientations, XfiQ 1 and ASwP*. 
We first consider [110] channels, used to determine A^fg 
[Figs. [2fa) and[2fc)]. For the long channel, with x id = 
25 /xm, the spin signal is suppressed below detectable 
levels between 4 and 6 T [Fig. [2|c)], indicating that 
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FIG. 2: Magnetic field dependence of V n i measured with in- 
jector and detector set to le 2 /h and separated by 5 /im (a, b) 
and 25 /im (c, d). Channel orientations are shown at the top: 
along [110] for (a, c) and along [1 10] for (b, d). In (c) and 
(d), dashed lines are V n i measured with A-gate undepleted 
[see Fig. [ljd)]. Finite signal at zero field in all panels is a 
signature of the Peltier effect, not of spin current '23 |. 



Aff R <C 25 //m but making it difficult to determine an 
accurate numerical value. (The shallow slope between 4 
and 6 T is a slowly increasing thermoelectric signal, and 
not related to spin [23].) For the short channel, with 
%id = 5 (im, the spin signal is detectable throughout the 
BSR dip, with a minimum at 6.5 T. V n i is similar before 
(B ext = 4 T) and after (B ext = 10 T) the dip in Fig. 2(a), 
indicating that injector and the detector QPCs are fully 
polarized for B ext > 4 T in this device. Thus changes 
in V n i ar e due only to changes in spin relaxation length 
caused by BSR, and A can be deduced from solutions to 
a ID spin diffusion equation, d 2 V n i/dx 2 — V n i/X 2 giving 
A n s R = 3.3i0.7^m[llj. 

Channels oriented along [110], on the other hand, 
showed a clear spin signal on resonance for both long and 
short geometries [Figs. 2(b) and 2(d)]. The signal did not 
saturate with field up to B ext =12T in these devices, pos- 
sibly due to QPC polarizations slowly increasing with 
field, so it was not possible to determine Xf^ simply 
by comparing V n i inside and outside the BSR dip. In- 
stead, the A-gate in the long channel was used to change 
the channel length in situ [Fig. 1(d)]: a polarization- 



independent strategy for extracting A BSR [23] • V n i de- 
creases when the A-gate is undepleted because the right- 
hand reservoir with equilibrium spin polarization is ef- 
fectively brought closer to the detector. Spin relaxation 
length is extracted from the ratio V n i (A-gate on)/V n i(A- 
gate off) using the solution to the ID diffusion equation, 
giving 65±5 fim off resonance and 40±5 fim in the BSR 
dip [24]. 




FIG. 3: (a) Monte Carlo simulations of spin relaxation time, 
converted into spin relaxation length A (see text), using pa- 
rameters I = 12 fim, vp = 8.4 x 10 4 m/s, and w = 900 
nm corresponding to the channel in Fig. 2(a). From top to 
bottom, curves represent (a — /3)=0.25, 0.5, 0.75, 1, 2, 3, 5, 
7, and 10 meVA. The value of (a + /3) is fixed at 1 meVA. 
(b) Spin relaxation length at the center of the dip, A BSR , for 
channel parameters in Fig. 2(a) (Auo R ) and Fig. 2(d) (Xf^) 
plotted against the relevant combination of a and /3. The re- 
lation ,\ BSR oc (a ± P)^ 1 valid in the weak spin-orbit limit is 
shown for reference (solid line). For the channel in Fig. 2(d), 
£ = 18 /im, vf = 9.9 x 10 4 m/s, and w = 950 nm. 

For Dyakonov-Perel' spin relaxation in which the spin 
precession length due to SOI is longer than the chan- 
nel width, a spin relaxation time of approximately r sr oc 
(B so )~ 2 is expected, corresponding to a spin relaxation 
length A = ^T sr £v F /2 cx (B 30 )" 1 where I is the mean 
free path (20] [25] . In order to obtain more accurate esti- 
mates of r sr , spin dynamics were simulated using Monte- 
Carlo techniques, averaging over an ensemble of realistic 
random trajectories based on the mean free path, £, chan- 
nel width, w, and vf for a particular channel [24 . The 
simulations provided estimates of T sr (B so ) for the two 
channels, which were converted into A and matched to 
the experimental data. 

An example is shown in Fig. 3(a), where Xno(Bf x ^) 
is obtained for a wide range of spin-orbit parameters 
(a — /3), from simulations based on t, Vp, and w for the 
channel in Fig. 2(a). The simulations confirmed that 
BSR in the 110 channel was sensitive only to (a — 0), 
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but nearly independent of (a + /?), and vice versa (c.f. 
Eq. [lj; the sum and difference of a and (3 could there- 
fore be fitted independently to the two channels, greatly 
facilitating the determination of both parameters. An 
analogous simulation was performed for the channel in 
Fig. 2(d). Based on this analysis, the values of A BSR 
found earlier gave (a-0) = 7±2 meVA and {a + p) = 
0.45±0.05 meVA: an anisotropy of \a-/3\/\a+p\ = 15±5. 
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FIG. 4: Ballistic spin resonance induced by out-of-plane mag- 
netic field B e ^ 1 in a [110] channel, (a), and a [110] channel, 
(b), both with Xid = 25 )im. Schematics on top show the 
channel and magnetic field orientations, and simplified cy- 
clotron trajectories for the two panels. In each case V„i in the 
dip must be compared to signal for B™' — 0, where no dip is 
seen, to take into account partial QPC polarization. 

BSR can also be induced with B ext oriented primar- 
ily along the channel axis if the component of B ext per- 
pendicular to the 2DEG plane, Bj^ f , is not zero [15] . 
Electrons follow cyclotron trajectories due to B ext , along 
which both ki w and fcno oscillate. This ensures an os- 
cillating component of B so perpendicular to B ext , as re- 
quired for BSR. Figure|4]demonstrates how BSR develops 
with B e ^ 1 in orthogonal channels. For B ext = the BSR 
dip is not observed because the only oscillating compo- 
nent of the momentum is perpendicular to the channel. 
The oscillating component of the spin-orbit field therefore 
lies along the channel axis, parallel to B ext . At higher 
B ext the BSR dip grows because the component of mo- 
mentum along the channel axis also begins to oscillate. 
The dip shifts to higher in-plane field with B ext : as the 
bouncing frequency increases with decreasing cyclotron 
radius. 

BSR from B e f f in [110] and [110] channels confirms 
the strong spin-orbit anisotropy reported above. The dip 
in the [110] channel is now weaker than in [110] channel, 
because Sffo probes the cyclotron-induced oscillation in 



fcno, the weak spin-orbit direction. In the [110] channel 
BSR leads to a suppression of V n i by nearly 85% com- 
pared to the signal for B ext = [Fig. 4(b)], while V n i is 
only suppressed by 35% in the [110] channel. 

In conclusion, BSR is used to estimate the degree of 
spin-orbit anisotropy in a GaAs/AlGaAs 2DEG. The spin 
relaxation length at BSR is found to be more than an 
order of magnitude shorter in channels oriented in the 
[110] direction compared to [110]. The extracted values 
for a and f3 are within 20% of each other, suggesting 
that in the future it may be possible to reach the regime 
|a| = \0\ by tuning the 2DEG with a top gate [15112^] . 
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